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ABSTRACT 

We analyze the density field of 264,283 galaxies observed by the Sloan Digital Sky Survey 
(SDSS)-III Bai-yon Oscillation Spectroscopic Survey (BOSS) and included in the SDSS data 
release nine (DR9). In total, the SDSS DR9 BOSS data includes spectroscopic redshifts for 
over 400,000 galaxies spread over a footprint of more than 3,000 deg^. We measure the power 
spectrum of these galaxies with redshifts 0.43 < z < 0.7 in order to constrain the amount of 
local non-Gaussianity, /^l^'^ in the primordial density field, paying particular attention to the 
impact of systematic uncertainties. The BOSS galaxy density field is systematically affected 
by the local stellar density and this influences the ability to accurately measure /nl^'- I" 
the absence of any correction, we find (erroneously) that the probability that /nl^' is greater 
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than zero, P{f'^ff' > 0), is 99.5%. After quantifying and correcting for the systematic bias 
and including the added uncertainty, we find — 92 < /^l^' < 398 at 95% confidence, and 
-P(/nl^' > 0) = 91.0%. A more conservative approach assumes that we have only learned 
the fc-dependence of the systematic bias and allows any amplitude for the systematic cor- 
rection; we find that the systematic effect is not fully degenerate with that of /nl**'- and we 
determine that -168 < < 364 (at 95% confidence) and P{fj^l''^ > 0) = 68%. This 

analysis demonstrates the importance of accounting for the impact of Galactic foregrounds 
on /^l'^' measurements. We outline the methods that account for these systematic biases and 
uncertainties. We expect our methods to yield robust constraints on /^l^' for both our own 
and future large-scale-structure investigations. 

Key words: cosmology: observations, large-scale structure 



1 INTRODUCTION 

Measuring the clustering of galaxies on large scales provides an ex- 
citing opportunity to test inflationary models. Slow roll, single field 
inflation is believed to generate a primordial gravitational potential 
(hereafter potential) that can be described statistically by a (nearly) 
Gaussian random field ( [Bardeen et al.|1986j . Alternative inflation- 
ary models (e.g., multiple field) predict there to be significant non- 
Gaussian components to the potential (see, e.g., |Waridsl|2010| for 
a review). For a Gaussian random field, the amplitude of fluctua- 
tions in 3-point configurations is zero. It is therefore convenient to 
express the degree of non-Gaussianity, /nl, as a dimensionless ra- 
tio between the amplitude of the bispectrum, Bi>{ki, k2, fca), and 
power spectrum, P$(fc), of the fluctuations of the total potentiaQ 



2 [P^(ki)P^{k2) + Pi,(k2)P.S.{kz) + P^(kl)P^{k3 



.(1) 



One can denote the portion of the potential that can be de- 
scribed as a Gaussian random field as and assume that /nl is a 
function of the potential (i.e., that it is 'local'). To 2nd order, this 
approach yields 



' + /Nr'(0'-(0')). 



(2) 



For this simplest treatment and evaluating Eq. [T|in the limit where 
fc-space triangle configurations satisfy lfci2l <Si \kri\, \ k23,\ (known 
as the "squeezed" limit), it has been shown that a bias in the dark 
matter halo power spectrum proportional to 



would result 



dPalal et al.l2008',' Matarrese & Verde|2008l l. [Komatsul ( [20Tol l and 
[Creminelli et al. ( 20TT| have shown that, for standard, single-field 
inflation, the amplitude of the bispectrum in the "squeezed" limit 
becomes vanishingly small, and thus any detected scale-dependent 
bias at large-scales would disprove all single-field models, subject 
to the condition that the field starts in the vacuum. Therefore, mea- 
surements of the large-scale clustering of galaxies have the poten- 
tial to yield significant insight into the physical mechanism that 
drove inflation. 

Primordial non-Gaussianity can also be measured from the 
bispectrum of cosmic microwave background (CMB) anisotropics 
(see, e.g., |Bartolo et al.|2004[|Komatsu|2010 and references therein) 
and, in principle, galaxi es (see e.g. |Sefusatti & Komatsu||2007[ 



aljg 



Scoccimarro et al. 20121. WMAP7 found 



10 < fW^< 74 to 



95% confidence i Komatsu et al.|2011^ . To date, the reported con- 



straints on /nl'' from galaxy power spectrum measurements have 

^ Specifically, Bardeen's gauge-invariant potential, which is equivalent to 
the Newtonian potential at sub-horizon scales 



been competitive with those obtained from the CMB bispectrum. 
[Slosar et al.| ( [2008] > used a combination galaxy and quasar cluster- 
ing measurements to find -29 < /nl''' < 70 at 95% confidence. 
|Xia et al!lj201 1| > analyzed similar measurements based on updated 
data to decrease the 95% confidence interval to 5 < /nl''' < 84. 
Komatsu et al.|p011| l combined their CMB results with the |Slosar| 

' < 59 at 95% confi- 
•^^l^'K but do not rule 
(2012 I predict that combin- 



et al. 



1 2008 I result to obtain -5 



^ ("local 

< /nl 

dence. Current r esults thus favour positive 



out /i°£"' 



0. 



Giannantonio et al. 



ing CMB data from the Planck mission ^The Planck Collaboration! 
|2006[ l and galaxy data from a Euclid-like ^Laureijs et al.|201 1 1 sur- 



vey will reduce the la uncertainty on /nl'' to 3. 

A large value of /nl''' implies larger amplitudes in the 2- 
point clustering of galaxies at large separations than expected in 
the standard ACDM paradig m. Recent studies (e.g., | Thomas et al.] 
|201 l[|Sawangwit et al.|201 l[|Nikoloudakis et al.|2012| l have indeed 
found larger than expected clustering amplitudes at large scales, us- 
ing photometric Sloan Digital Sky Survey (SDSS; [York et al.|2000l > 
data. However, iRoss et al. , 1 ,201 Ij showed that the excess found in 
[Thomas et al.| (|20TTJ was due, at least partially, to stellar contam- 
ination, and that systematics correlated with the Galaxy (e.g., stel- 
lar density and Galactic extinction) may impart spurious clustering 
signal at large scales. [Nikoloudakis et aT] l |2012[ l find their mea- 
surements yield /nl"' = 90 ± 30 at 68% confidence, but suggest 
this result may be better interpreted as /nl"' < 120 at 84% confi- 
dence, in light of the systematic concerns. General relativistic (GR) 
corrections are also expected to cause features in the clustering of 
galaxies at the largest (super-horizon) scales. These effe cts are ex 
pected to be small compared to that of /nl^'^ as, e.g. 



et al. 



Maartens 



1 2012 1 find that the effects of /nl" on the power spectrum 



dominate GR corrections for /nl"' > 10. 



We analyze the SDSS data release nine (DR9;Ahn etal."2012[ 
Baryon Oscillation Spectroscopic Survey (BOSS; ^Dawson et al 
|2012| > "CMASS" sample of galaxies. This sample has the largest ef- 
fective volume fTegm ark & Peebles|1998| l of any spectroscopic sur- 
vey and comprises approximately 1/3 of the final (planned) BOSS 
CMASS sample. The clustering of this data set has been extensively 
studied ( [Anderson et al.|2012|[Nuza et al.|2012[|Reid et al.|2012[ 
[Sanchez et al.[[20T2 Tojeiro et al.[[2012) . In particular, the sam- 
ple has been simulated with 600 mock galaxy catalogs ( [Manera et] 
[al.|2012| l and the issues identified by Ross et al. ( 2011) have been 
addressed via the application of an un-biased weighting scheme 
( [Ross et al.|2012f . We thus have the tools to robustly investigate 
the information content in the large scale clustering, even in light 
of systematic concerns. 
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Despite having the largest effective volume of any spectro- 
scopic study, the DR9 CMASS sample is smaller than existing pho- 
tometric redshift samples of galaxies and quasars (existing quasar 
samples have roughly 20 times the volume). Given that /nl^' im- 
parts a change in the expected clustering measurements that is most 
pronounced at large scales and is smoothly varying, we should not 
expect /nl^' constraints from the DR9 CMASS sample to be com- 
petitive even with those one could obtain with the SDSS data re- 
lease eight I Aihara et al.|201 1\ photometric redshift sample created 



in |Ross et a 



1 201 1 1 (which has approximately 3 times the angular 



footprint of our data set and was trained using early BOSS red- 
shifts), let alone a photometric quasar sample. However, we expect 
our study to have best-quantified the impact of systematics, given 
the analyses of |Ross et al.| ( |2012j l. 

We primarily use measurements of the spherically-averaged 
power spectrum, P{k), to constrain /nl''', as the scales most af- 
fected by local non-Gaussianity are most isolated in fc-space. In 
Appendix [C] we test for consistency using the spherically aver- 
aged correlation function, £,{s). The main purpose of this study is 
to demonstrate how the ability to constrain /nl^' using the galaxy 
power spectrum depends on the treatment of systematic uncertain- 
ties related to stellar density. We describe the observed and simu- 
lated data we use in Section|2]and our analysis methods in Section 
[5] Our results, and their dependence on our treatment of potential 
systematics, are presented in Section|4] In SectionBJwe discuss our 
results in the context of current and future /nl^ measurements, 
and we summarize our conclusions in Section|6] 

Unless otherwise noted, we assume a flat cosmology with 
= 0.285,^6 = 0.0459, /i = 0.70, = 0.96, and as = 0-8, 
as are approximately the best-fit values found by [Sanchez et al.| 
1 20T2I when fitting the full shape of ^(s). 



2 DATA 



The SDSS-III Baryon Oscillation Spectroscopic Survey ( |Eisenstem| 
|et al.|20n^ obtains targets using SDSS imaging data. In combina- 
tion, the SDSS-I, SDSS-II, and SDSS-III surveys obtained wide- 
field CCD p hotometry (G unn et al. 1998! [2006) in five passbands 
(u, g, r, i, z; |Fukugita et al.||I996) , amassing a total footprint of 
14,555 deg'^, internally calibrated using the 'uber-calibration' pro- 
cess described in [Padmanabhan et al.[p008] >, and with a 50% com- 
pleteness limit of point sources at r = 22.5 ( [Aihara et al.|20Il] l. 
From this imaging data, BOSS has targeted 1.5 million massive 
galaxies, 150,000 quasars, and over 75,000 ancillary ta rgets for 



spectroscopic observation over an area of 10,000 deg^ dDawson 



[et al.[[2012| l. BOSS observations began in fall 2009, and the last 
spectra of targeted galaxies will be acquired in 2014. The BOSS 
spectrographs (R = 1300-3000) are fed by 1000 optical fibres in a 
single pointing, each with a 2" aperture i Smee et al.|2012| . Each 
observation is performed in a series of 15-minute exposures and 
integrated until a fiducial minimum signal-to-noise ratio, chosen to 
ensure a high redshift success rate, is reached. Redshifts are deter- 
mined as described in |Bolton et al.| ((2012 1. 

We use data from the SDSS-III DR9 BOSS CMASS sample 
of galaxies, as defined by |Eisenstein et al^p011| l. We use the same 
DR9 CMASS sample and treat it in the exact same way as in "An-' 
[derson et al.lpOTI) and [Ross et aLlpOnj l. This sample has 264,283 
galaxies spread over an effective area of 3275 deg^ 2584 deg^ of 
which is in the North Galactic Cap. 

We use the 600 mock DR9 CMASS catalogs (henceforth 
"mocks") generated by [Manera et al.[ ( [20I2[ l to generate covari- 



ance matrices for clustering estimators. [Manera et al.[ |20T2| used 
the initial conditions of a flat cosmology defined by Sim = 
0.274, Slft/i^ = 0.0224, /i = 0.70, = 0.95, and as = 0.8 



(same as the fiducial cosmologies assumed in [White et al.[[20Il[ 
and |Anderson et al.|20I2| l and generated dark matter halo fields at 
redshift 0.55. These simulations were produced using a 2nd-order 
Lagrangian Perturbation Theory (2LPT) approach inspired by the 
Perturbation Theory Halos (PTHalos; Scoccimarro & Sheth 2002]l 
model. Galaxies were placed in halos using the halo occupation 
distribution determined from measurements of the correlation func- 
tion of CMASS galaxies and the parameterization of [Zheng et al.[ 
( [2007) 1. The DR9 angular footprint was then applied and galaxies 
were sampled along the radial direction such that the mean n{z) 
matched the CMASS n(z), thereby providing 600 catalogs simu- 
lating the observed DR9 CMASS sample. See [Manera et al.[j2012l > 
for further details. The mocks have /nl^' — and we describe how 
we account for this in our covariance matrices in Sectionl3^ 



3 ANALYSIS TECHNIQUES 



3.1 Physical Model 



As shown in, e.g., [Matarrese et al.[ j2000[ l, primordial non- 
Gaussianity alters the expected mass function of dark matter halos, 
with positive /nl^' yielding more high-mass dark matter halos. As 
outlined in both Dalai et al. ( 2008) and [Matarrese & Verde[ ( |2008| l, 
this causes a scale dependent bias in the power spectrum of dark 
matter halos. These studies have shown that, expressing the Gaus- 
sian halo bias as b, the additive bias is given by 



where 



(3) 



(4) 



and Sc{z) = 1.686/ D{z) is the critical spherical over-density for 
dark matter halos to collapse at redshift z, D{z) is the linear growth 
factor (normalized to equal 1 at 2 = 0), and T(k) is the transfer 
function. We obtain the z = linear power-spectrum, Puiik), and 
transfer function using the software package CAMEPj ^Lewis et al. 
[2000| l and our fiducial cosmology. 
The total bias is given by 



btot = & + BML{k). 



(5) 



Thus, the galaxy power spectrum, Pg{k) can be expressed given 
the matter power spectrum, PM{k), as 

Pg{k) = b^Puik) + 2bBNL{k)PM{k) 

+BlL{k)PM{k). 



(6) 



We will measure spherically averaged clustering in redshift space. 
Assuming linear theory, the redshift space, galaxy power spectrum, 
Pg{k) is determined from the real space galaxy power spectrum 
( [Kaiser|1987") 



P°{k) = Pg{k) + 2/3fbtotPM{k) + l/5,fPM{k), 



(7) 



where / = dlnD/dlna. The redshift zero result is then evolved 
via Pg{k, z) = D^{z)Pg{k, 0). We restrict our analysis to fc < 



see the website camb.info 
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0.05/iMpc~^, where our linear model is expected to be a valid ap- 
proximation, as is justified via comparison to the results of mock 
DR9 CMASS samples described in Section [7!2l 

Various studies (e.g., Desjacques et al. 2009{|Wagner & Verde| 



2012 Scoccimarro et al. 2012 1 have shown that non-zero /^l 



produces corrections to the scale-dependent and -independent halo 
bias that are more complicated functions of halo mass (and age) 
than the simple model we have presented above. Further, Des^ 
[jacques et aL]^2011^ and |Roth & Porciani| ( (20T2 1 show that higher- 
order contributions to the (non-Gaussian) primordial potential may 
bias determinations of /nl^'. These concerns suggest that measure- 
ments of /nl^' using the techniques we describe above may be 
biased and that any precise determination of /nl''' must consider 
these effects. However, in all of the models, detection of non-zero 
/nl''' should occur only when there is local non-Gaussianity in the 
primordial potential. The values we determine should be viewed as 
effective /nl^' values, but we expect the relative importance of sys- 
tematic biases imparted by observational effects to be independent 
of these model concerns. 



3.2 Measuring the Power Spectrum 

We measure the power spectrum, Pm{k), using the standard 
Fourier technique of Feldman et al. ( 1994), as described in |Reid| 
|et al.|(2010) and ,Anderson et al., ( ^2012^ . In particular, we calcu- 
late the spherically-averaged power in k bands of width AA; — 
0.004/!,Mpc~^ using a 2048"^ grid. For all measurements we 
present, the galaxies are weighted to account for missing close- 
pairs (the BOSS spectrographs cannot simultaneously observe ob- 
jects separated by less than 62") and redshift failures. We also ap- 
ply weights to account for the systematic relationship between tar- 
get galaxy density and stellar density, and compare to results when 
these weights are not applied. We label measurements made using 
the stellar density weights as Pm,star, and treat these as the fidu- 
cial measurements. We label those made without the stellar density 
weights as Pm,nw Finally, we apply 'FKP' (standing for Feldman- 
Kaiser-Peacock) weights to both galaxies and randoms using the 
prescription of |Feldman et al.|{i994^ , which amounts to a redshift- 
dependent weighting in our application. The total weights for each 
galaxy/random are summed at each grid point. The process of cal- 
culating weights is described in detail in |Ross et al.| ( |2012| l. 

In the resulting measurements, k bands will be correlated as a 
result of the survey geometry. We measure this effect by determin- 
ing the power spectrum of the survey geometry, as simulated using 
an unclustered random catalog. This exercise yields the spherically- 
averaged "window power", which we denote Pwin{k). To com- 
bine the windows of the NGC and SGC, we determine the volume- 
weighted average of the respective Pmin{k). 



3.3 Accounting for the Survey Window 

The measured Pm(k) is a convolution between the window power 
and the true, underlying, power spectrum, Pt{k). In practice, we 
measure Pm{k) in discrete A; -bands. Thus, we determine a "win- 
dow matrix", Ty[fci][fcj], using the spherically-averaged Pwin and 
the process outlined in Appendix [A| Our expectation for the mea- 
sured Pconv (k) is given by 




k (hMpc-i) 



Figure 1. The mean power spectrum recovered from the mocks, Pm(k), 
divided by the input power spectrum used to generate the mock galaxy cat- 
alogs that has been first translated to redshift space at z = 0.55 assuming 
a real-space bias of 1.889 (the best-fit value) and then convolved with the 
window applied to the mock footprint. The error-bars are the 1 cr uncertain- 
ties on the mean of the 600 mocks and the dashed lines reflect the standard 
deviation of the 600 mocks. 



(9) 



(8) 



The second term, which we will here after refer to as the "window 
subtraction", in Eq.jsjis necessary because the measured Pm{k) is 
zero at fc = by definitiorj^ Inspection of Eqs.js] and [9] reveals 
that the normalization given by Po yields Pconv (0) ~ 0, and thus 
Pconv{0) ~ Pm{0) = by construction. We measure Pmin and 
calculate Pg{k) in bins AA: — 0.0005/iMpc~^, yielding sufficient 
resolution to calculate Pconv (k) in (the measured) bin width AA — 
0.004/iMpc-^ 

We test our determinations of Pwin and W^[A;i][Aj] using 
Pm(k) measured from 600 mock galaxy catalogs generated by 
|Manera et al.|pOT2) . We should be able to recover the mean of 
the Pni{k) measured for the mocks by convolving the input P{k) 
used to create the mocks and the P^in measured from the win- 
dow applied in |Manera et al.|f2012| >. The ratio between the mean 
Pm(k) recovered from the mocks and the convolved input power 
spectrum, Pconv{k) (translated to redshift space at 2 = 0.55 and 
assuming the cosmology used to generate the mocks and a real- 
space bias of 1.891), is displayed in Fig.[T| The error-bars reflect 
the uncertainty in the mean, and the dashed lines display the stan- 
dard deviation (which is the expected statistical uncertainty for our 
CMASS measurements). We determine a best-fit bias of 1.891 by 
scaling the full covariance matrix (i.e., we divide each element by 
600). Only for the bin centred at A = 0.002/iMpc"^ does the re- 
sult differ by more than la from this best-fit theoretical expecta- 
tion. However, the correlation between fc-bins is significant and the 
minimum is 21.1 (11 degrees of freedom and 12 fc-bins). The 



where 



This is analogous to the integral constraint on the correlation function. 
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Figure 2. Top panel: The input redshift space power spectram at z = 0.55 

, di- 



vided by P°(fc) for /ij 



0. Bottom panel: The same information as 



the top panel, except the models have now been convolved with the window 
function and the range of the P°(k)/ P°{k, /^l^' = 0) axis has been 
decreased by more than an order of magnitude. 



value is small enough to suggest our measurements of /nl^' should 
be unbiased. Indeed, we find f^l""^ = 1.3 ± 5.0 when we fit the 
average of the mock P{k). Our treatment of the window therefore 
introduces no measurable bias on the recovered /nl^'. This uncer- 
tainty of 5.0 corresponds to the mean of the 600 realizations of the 
DR9 CMASS sample; we therefore should expect a la uncertainty 
on f^i;^^ close to 120 (~ 5 x ^600) for the actual DR9 CMASS 
sample. 

The window, and the associated window subtraction term in 
Eq.[9] imply that the expected measurement of the power spectrum 
will not have the same shape as the true underlying power spectrum. 
This effect is illustrated by Fig.|2] In the top panel, a series of input 
Pg{k) (evaluated at the centre of the Afc = 0.0005/iMpc"^ bins) 
are plotted for -300 < /nT"' 300, divided by the input P°{k) 
with = 0. At the lowest k, negative /^f''' actually causes 



an enhancement of the power spectrum; this is due to the fact that 
the term proportional to the square of /nl^' dominates and thus the 
square of the total bias becomes larger than in the case of /nl''' ~ 



0. The bottom panel displays the same information as the top panel, 
after Pg{k) has been convolved with the window (now for Ak = 
0.002 bins). At the low k bins, the effect of /nl'"' is decreased 
by more than an order of magnitude when compared to the input 
Pg{k). At these values of k, the measurement depends strongly 
on the underlying Pg{k) at larger k and Pmin, and their ability to 
constrain /nl''' is much weaker than one naively expects. 



3.4 Modelling the Systematic Contribution 

We compare the measured power spectrum, Pm{k), of CMASS 
galaxies to model Pm.od{k, /nl'^'i Si), which combine the theoret- 
ical galaxy power spectrum convolved with the CMASS window 
function and a term 5* (described below) to model the contribution 
of systematics to the measured power. |Ross et al.H2012^ found that 
a weight for stellar density had a large effect on the measured low-fc 
CMASS power spectrum. We assume that the correction for poten- 
tial systematics always has the same form, but the magnitude of the 
correction can vary. This assumption is justified by the fact that we 
found this to be true when testing different weighting schemes in 
|Ross et al.| ( |2012) . Our total model is 



/J /"local o\ 
iod(K, /nL j'J) 



Pconvik, f'H'L )~\~S\Pjri 



■j{k) Pm,star{k)^^ 



where Pm.starik) denotes the measurement made using the stel- 
lar density weight, and Pm,nw{k) denotes the measurement deter- 
mined when the stellar density weights are not applied. Treating 
Pm,star as the fiduclal measurement, this approach fixes the sys- 
tematic correction in the observed power spectrum and uses the 
model to account for the possibility that our measurement has a 
remaining systematic bias by scaling the amplitude of the fiducial 
correction. Specifically, S — represents the fiducial case assum- 
ing we have properly removed any systematic bias from our mea- 
surement and positive/negative S accounts for the possibility a sys- 
tematic bias remains/has been incorrectly subtracted by our treat- 
ment of the data. 

We apply weights to account for the systematic effect of stars 
when we determine the power spectrum. These weights were deter- 
mined via a linear fit to the relationship between the target galaxy 
density and the stellar density. We can estimate the uncertainty on 
5* by determining these weights for each of the 600 mocks (which 
have no intrinsic systematic relationship). This analysis was per- 
formed in Ross et al. 1 20121, where ^(s) was measured for each 
mock sample with weights determined and applied in each individ- 
ual case. The weights did not bias the mean of the measurements or 
the variance (see their Fig. A2), but each individual measurement 
was given a nearly constant bias with standard deviation 1.5x 10~*. 
This standard deviation on the bias is 10% the size of difference 
between the unweighted, £,„w, and weighted CMASS £,{s). There- 
fore, we expect the statistical uncertainty on 5 to be 0. 1 for the DR9 
CMASS sample, and we believe allowing any value of 5* accounts 
for additional, unknown systematic effects. 



3.5 Comparing Measurement to Model 

Previous DR9 studies (e.g., 'An derson et al.|2012 t'Ross et a l.|2012| l 
used measurements of the power spectra of 600 mock realizations 
of the CMASS sample to construct the covariance matrix. These 
mock realizations are for /nl^' — 0. For a constant number den- 
sity, n, the total variance is expected to scale as (see, e.g., |Feldman| 
|etal.|1994| i 



(jp oc {P+l/nf. 



(11) 
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The value of /nl^ significantly alters the expected power spec- 
trum at low k, and thus we expect the amplitude of the covariance 
to depend significantly on the true value of /nl^'. The number den- 
sity for the DR9 CMASS sample is ~ 3 x 10"*/i^Mpc"^ and its 
measured power spectrum is greater than 2.5 x 10*ft~''Mpc'' across 
all of the scales we consider Thus, we expect cosmic variance to 
dominate the uncertainty at these scales by a factor of at least 8 
(the ratio of P to 1 /n). Therefore, the expected uncertainty on the 
power spectrum of can be approximated by crp oc P and thus the 
expected uncertainty on a^/P can be approximated as independent 
of the theoretical model. 

The uncertainty on ln[P] is Therefore, and based on the 
arguments above, we compare ln[Pm] to ln[P,„od] and thus when 
testing models against our measurements, we minimize the 
given by 



Pmod (^i ) 



Pmod {J^j ) 



(12) 



where Cp ^ is the covariance matrix of ln[P,n]. Cp ^ is determined 
from the mocks via 



Cj> \ki\ \kj 



^ 600 



599 



Pm iki 



In 



Pm{kj) 



(13) 



where P," represents the measured power spectrum of mock n and 
Pm represents the mean across all mocks. 

For the linear theory power spectrum, the correlation between 
different fc-bins is the effect of the survey window alone. Thus, 
in the cosmic variance limit (in which the diagonal elements are 
constant), the covariance matrix of the logarithm of the power for 
k < 0.05/iMpc~^ depends primarily on the survey window. Thus, 
to the extent that we can ignore model dependent variations in the 
shot-noise term and non-linear effects on the covariance of ln[P] 
for k < 0.05/iMpc^^, we expect our treatment to be correct in- 
dependent of /nl^\ or any other cosmological model that changes 
the expected amplitude of the power spectrum on large scaleQ Our 
choice to test ln[P(fc)] is further supported by the fact that, as de- 
scribed in Appendix[B| we find that the skewness of the distribution 
of ln[P(fc)] is significantly smaller than that of P(k). 



We are using the same Pm (fc) measurements as Anderson et 
[ar] ( |2012} , which assumed a flat geometry with Sim = 0.274, but 
we adopt a flat geometry with Q,m = 0.285 as our fiducial cos- 
mology. In order to account for the difference, we assume spher- 
ical symmetry and dilate the center of the fc-bin by D„(S7m = 
0.285)/D.„(S7m = 0.274) = 1.006, where 



= \cz{l 



zfDlH 



-ni/3 



(14) 



where Da is the angular diameter distance and H is Hubble pa- 
rameter. Thus, when comparing the theoretical power spectrum 
to our measurements, we determine the input Pg{k) in fc-bins 
of width 5.03xlO"'*/iMpc~^ and the convolved Pconv{k) in fc- 
bins of width 4.024X lO^'^/iMpc^^ (whereas the fc-bin sizes were 
5xlO"*/iMpc-^ and 4xlO"^ftMpc"^ for the flat ilm = 0.274 
cosmology). Dilations are applied in the same manner when alter- 
native values of Qm are tested. 



* Our covariance matrix does include both shot-noise and non-linear ef- 
fects, it simply will not account for the possibility that these effects scale 
differently with changes in, e.g., than the power 



0) 0.04 
N 



'(i) fiducial (P„,,,„) ■ ■ ' ■ ■ 

(ii) no syslematics correction (P^ 

(iii) Gaussian prior S = 0±0.1 

(iv) marginalizing over S(P„„,-P„ 




f local 

Figure 3. The noiTnalized (so that they integrate to 1) probability distri- 
butions for the local non-Gaussianity parameter /^l^'' treat- 
ment of systematics, applied to the DR9 CMASS sample. The blue curve 
shows the result using our fiducial treatment which uses the power spec- 
tram, Pm.star determined using the stellar density weights; the green 
curve shows the result when we use the power spectram measurement, 
Pm,nw, that does not include the stellar density weights; the light blue 
curve shows the result when we marginalize over an additional term 
S{Pm,nw — Pm,star) in the model and allow it vary within a Gaussian 
prior of it 0.1; and the red curve shows the results when we marginalize 
over S and allow it vary freely. 



4 RESULTS 

We present results for four treatments of the systematic contribution 
to the model described by Eq.|10| 

(i) S = 0; this corresponds to our best estimate of the power 
spectrum, measured using the weights for stellar density, but we 
ignore systematic uncertainty in the application of these weights; 

(ii) S = — 1; this corresponds to the case where no weights are 
applied for stellar density; 

(iii) S is allowed to vary around with a Gaussian prior of width 
0.1 (this value is as motivated by the tests performed on mocks in 
'Ross et al. 2012) accounting for the uncertainty in application of 
the stellar density weight; and 

(iv) Any value S > — 1 is permitted. This final case assumes we 
have learned the fc-dependence of systematic effects on our mea- 
surement, but allows for the possibility of an unknown systematic 
(S > 0) or that the stellar density weights were an over-correction 
(5 < 0). 



We marginalize over S for cases (iii) and (iv) and b in all cases. 
There is a small degeneracy between Qm and /Nf'. We find that 
the best-fit value of /nl^' shifts by 12 comparing case (i) with 
Qm ~ 0.26 and Qm ~ 0.3. This change is approximately 10% 
of the expected statistical uncertainty, and is thus small enough 
that we expect the conclusions of this study are robust to changes 
in the background ACDM cosmology (flat ACDM with Qm = 
0.285, Qt = 0.0459, h = 0.70, = 0.96, and erg = 0.8). 
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Table 1. The recovered fj^{^ constraints and quality of fit for the four treatments of the systematic contribution to the model that we consider, where C.l. 
stands for confidence interval, ML stands for maximum likelihood, and < 0) is the sum of the probability at /^^^^ less than zero. 



Case 


jk,cal 95,7^ J 


If^l^^l 95% C.I. 


f^r' ML 


overall best-fit f^l'^'- 


xLst-fit/dof 


•S'bcst-fit 




i 


(-55, +400) 


(0, 375) 


217 


215 


15.9/10 


= 


5.9% 


ii 


(+65, +405) 


(65, 405) 


250 


250 


32.4/10 


= -1 


0.5% 


iii 


(-92, +398) 


(0, 365) 


209 


215 


15.6/9 


0.08 


9.0% 


iv 


(-168, +364) 


(0,315) 


155 


-98 


13.0/9 


0.45 


32% 




k (hMpc-') 



local = and 5 = is displayed 
Pm.star) from the measured 



Figure 4. Top panel: The measured DR9 CMASS P{k) for the labeled 
treatment of systematics (points with error-bars) and the associated best-fit 
model (solid lines). The model with fj^-^ 
with a dotted line. We subtract S{Pm,nw 
power spectrum, rather than add it to the theoretical model; in this way 
the points show the measurement assuming the given systematic treatment 
is the conect one. Bottom panel: The difference between the logarithms 
of a given power spectium and the overall best-fit model power spectrum 
(which is for S = 0.45, f^£^^ = —98), where the black circles represent 
the measured power spectrum (using weights for stellar density, Pm.star) 
and the lines represent the same four models and use the same scheme as in 
the top panel. 



Including the weights for stellar density, case (i), reduces the 
minimum to 15.9 (10 degrees of freedom), and we expect a 
value at least as large for 10% of consistent samples. The 
Axbcst-flt = 16.5 between cases (ii) and (i) shows that the stel- 
lar density weights are strongly prefeiTed. However, the maximum 
likelihood value of /n^' decreases only from 250 to 215. On can 
see in the top panel of Fig.|4]that the convolved models using these 
two /nl''' values (plotted so that the systematic term is applied to 
the measurement, not the model) appear qualitatively similar. The 
measurements at scales k > 0.02/iMpc~^ (where the stellar den- 
sity weights are not important) do not accommodate a significantly 
larger it is clearly the tension between the power spectrum 

measurements in k bins greater and less than k — 0.02/iMpc^^ 
that yields such a large in case (ii). 

The comparison between cases (i) and (iii) illustrates the im- 
portance of including the uncertainty on the systematic correc- 
tion. In both cases, we use the measured power spectrum that in- 
cludes the weights for stellar density. In case (iii), we are essen- 
tially adding a systematic uncertainty, by allowing S to vary around 
with a Gaussian prior of width 0. 1 (which is the uncertainty we 
determine from DR9 mocks), while for case (i) S is simply set to 
zero, meaning we have applied a systematic correction but allowed 
no uncertainty due to the correction. Comparing the PDF of case 
(iii) to that of case (i), we find that the PDF widens such that the 
95% confidence interval increases by 35 (representing an 8% in- 
crease), its peak shifts by A/^l''' = 8, and the probability that 
/nl'"' < increases from 5.9% to 9.0%. While these changes are 
not dramatic, they are significant and suggest (assuming the uncer- 
tainties add in quadrature so that a^y^ — at^t ^ f^atat) that the 



1(7 systematic uncertainty on /nl^ 
~45. 



for the DR9 CMASS sample is 



Fig.[3]displays the recovered probability distribution functions 
(PDFs) when we fit our power-spectrum measurements using the 
four treatments of the systematics we consider, and Table[T]lists the 
relevant attributes in each case. Only for case (ii), where we use the 
measurement without the weights for stellar density, (displayed in 
green) does the probability distribution appear Gaussian. For this 
case, /nl^' < is allowed at 0.5% and the PDF peaks at 7^;^°,^^' = 
250. In the absence of any systematic correction, we would have 
appeared to detect non-zero /nl^'. However, the value of the at 
the minimum, Xbcst-flt = 32.4 (10 degrees of freedom), for case 
(ii) indicates a problem as only 0.03% of samples consistent with 
our model would yield as large a x^- This result, in and of itself, 
suggests the analysis performed in case (ii) has a systematic in the 
treatment of the data/modelling. 



In case (iv), we allow any value of S. As shown in Fig.|3] 32% 
of the PDF (displayed in red) is now at < 0. Further, there 

are two peaks; the most likely value is at /nl''' = 155, but there is 
a second peak at f^l'^^ = -64 with an amphtude that is 79% that 
of the overall peak. As described in Appendix [B] when case (iv) 
was tested on the 600 mock realizations, multiple peaks are found 
in 24% of the PDFs (and in only 9% of the PDFs for the other 
cases). The minimum x^ for case (iv), 13.0 (9 degrees of freedom), 
is the overall minimum across all cases. The minimum occurs when 
5" = 0.45 and /Jl?/;^' = -98. Inspection of Fig. [4] reveals that the 
negative /nl^' value allows a better fit to the two measurements 
over 0.024 < k < 0.03ftMpc^^ and the 45% stronger correction 
allows the measurements at large-scales to be in good agreement. 
The x^ is reduced by 2.9 with the addition of one extra parame- 
ter (comparing to case i), and we would expect 16% of consistent 
samples to have a larger x^- 
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5 DISCUSSION 



5.2 Lessons for the Future 



5.1 The Meaning of Our Results 



The exact treatment of systematics causes large differences in the 
allowed range of /nl^', as the width of the 95% confidence in- 
terval increases from 340 to 532 from the least (case ii) to most 
(case iv) conservative treatment of systematics. The results of |Ross| 
|et al.|p012^ suggest a 10% uncertainty on the magnitude of the 
systematic correction applied for stellar density, and our case (iii) 
takes this effect into account, representing our treatment of system- 
atics including everything we know. The 95% confidence interval 
is -92 < < 398. 

A more pessimistic approach assumes that we know only the 
way that the systematic correction scales with k, but not its ampli- 
tude. This approach is treated in case (iv) and allows for the pres- 
ence of an unknown systematic. The 95% confidence interval for 
this situation is -168< /nl''' <364. The width of the 95% confi- 
dence interval increases by only 9% with respect to case iii. This 
result implies that more extreme values of /nl'^' are rejected by 
the P(fc) measurements at the scales relatively unaffected by the 
systematic corrections. The tests of our methods on the mocks (in 
Appendix [B} suggest that, for the DR9 CMASS sample, the confi- 
dence interval of the absolute value of /nl^' is most robust when 
less than 90% of the total probability is contained on one side of 
zero (i.e., 10% < P{f^t' < 0) < 90%). We find i/Jf^"'! < 315 
at 95% confidence for case (iv). 

The minimum x'^ decreases by 2.9 when we allow S to freely 
vary (the minimum is S = 0.45) compared to when we set S = 0, 
suggesting that there there may be a systematic in the treatment 
of the data that is undiscovered. We do not find that any smoothly 
varying model for P° can produce a similar decrease in the x^- 
For example, if the index of the scale-dependent bias is allowed to 
vary freely, we find the minimum x^ decreases by only 1.0 (when 
the scaling is A6 oc fc~^'*). We have taken a closer look at all of 
the potential systematics studied in |Ross et al.| l |2012) and found no 
further systematic effects on the measured clustering. We have also 
investigated the evolution of the clustering amplitude in the sample 
(a strongly evolving amplitude would preferentially increase the 
clustering at large scales) and found the variation to be less than 
10%, implying it has a negligible impact on our measurements. We 
suggest that this issue be revisited in future BOSS data releases, 
where the preference for a stronger correction should become more 
significant if there is an undiscovered systematic and there will be 
more data with which to test against potential systematic effects. 

Other studies have obtained much tighter constraints on /nl"' 
using measurements of galaxy/quasar clustering: JSlosar et al.| 



found -29 < 



rlocal 

/nl 



< 70 and 



5 < /nT'"' < 84 (both at 95 



Xia et al. 



(2011b found 



confidence). Both studies include 
SDSS photometric luminous red galaxy (LRG) samples that have 
a similar redshift range to that of our SDSS DR9 CMASS sam- 
ple (and significant overlap in angular footprints). Using only the 
SD SS ph otometric LRG data sample defined by |Padmanabhan et] 

268 < /J?L^' < 164 while 
NL^' < 351 using the LRG 



al 



j2007|l 



Xia et al. 



Slosar et al. 



. ( |2008| l found 

1 2011} measured -81 < /] 

sample created in [Thomas et al.| ( |2011^ . These values are similar 
to our own results and it is therefore doubtful that substituting our 
own for them would yield a significant change in the analyses that 
combine all samples. The final BOSS sample will be approximately 
three times as large as the DR9 sample and thus nearly halve (l/\/3) 
the statistical uncertainty on /nl^' that we find. 



On its own, the BOSS DR9 CMASS sample does not allow signif- 
icant improvement on existing /nl''' constraints. However, future 
studies can emulate our treatment of systematics uncertainties in 
order to obtain robust results. Here, we find that simply account- 
ing for the 10% uncertainty (implied by [Ross et al.| |2012 l in the 
correction we apply for stellar density increases the width of the 
95% confidence interval we determine for /Nf' by 35, implying 
fays = 45. It is not clear how the systematic uncertainty will scale 
in the presence of either more BOSS data or other data sets, e.g., 
this will depend on the exact distribution of stars across the angu- 
lar footprint. It is therefore possible that the systematic uncertainty 
will become, relatively, more (or less) important for future CMASS 
data samples or alternative large-scale structure probes. We recom- 
mend future studies subject their data samples to the level scrutiny 
presented in [Ross et aT] ( |2011[ |20I2^ . We expect any systematic 
dependency on a Galactic foreground (e.g., stars. Galactic extinc- 
tion, synchrotron emission) to impart a systematic bias and uncer- 
tainty onto /nl^' measurements; we caution that all future studies 
of large-scale clustering should test against this possibility in order 
to account for potential systematic errors. 

Concerns regarding systematics are especially relevant to the 
photometric quasar data that, to date, represent the large-scale 
structure da ta sample that yields the best /iSjL^' constraints, as 



etal. 



Xia 



( 201 1 1 find a 1(7 uncertainty on /nl'' of 26 when using only a 



SDSS data release six (DR6; |Adelman-McCarthy et al.pOOS) pho- 
tometric quasar sample. Photometric quasar samples are ideal for 
/nl^' studies, as the quasars occupy a large volume (decreasing the 
statistical uncertainty on large-scale measurements), have a large 
bias (increasing the relative effect of /nl'''), and are at high red- 
shift (and therefore have a smaller growth factor, which also in- 
creases the relative effect of /nl'''). They are also known to be 
subject to stellar contamination: [Myers et al.]j2007) found a stellar 
contamination of 4:A'^W%, for a photometric SDSS data release 
four (lAdelman-McCarthy et al."2006^ quasar sample constructed 
using the same methodology (Richards et al. 2004) as the quasar 
samples used by both Slosar et al. ( 2008 ) and Xia et al.| (201l) . 

The la bounds, using only quasars, were found to beSZ^^ 
[Slosar et al.lpOOSl l (for their fiducial QSO case using SDSS data 
release three (DR3; [Abazajian et al.|2003^ quasars) and 62 ± 26 
in |Xia et al.H201l| (using SDSS DR6 quasars). While the f^^:"-^ 
constraints are impressive, the two results are discrepant at ~ 1.5a 
(while using much of the same data). [Slosar et al.] ( |2008^ rejected 
quasars with photometric redshifts z < 1.45, after finding the cross 
correlation of this sample with stars had a non-zero amplitude. |Xia| 
|et al.| ( (20TT ) made no such cut on their quasar sample since they 
found no evidence of different contamination when splitting the 
data at z — 1.45, and they suggested that the difference was due 
to better calibration of the DR6 data compared to the DR3 sample. 
The degree to which stellar contamination could explain the tension 
between the two results is studied further in Giannantonio et al. (in 
prep.), where methods similar to our own are used to fully account 
for the systematic uncertainty of the effect of Galactic foregrounds 
on the measured clustering of quasars (and other tracers) and to 
obtain robust constraints on /nl"'- 



6 CONCLUSIONS 

• We have described a method that quantifies both the systematic 
bias imparted by Galactic foregrounds on clustering measurements 
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and its uncertainty, thus allowing one to obtain un-biased /nl'' 
constraints with realistic error estimates. 

• We find no reliable evidence for non-zero /nl'''. 

• The data show an extreme preference for our fiducial system- 
atic correction (the application of weights for stellar density when 
calculating the power spectrum, case [i]) compared to the no sys- 
tematics correction case (ii): the difference in the minimum is 
16.5, when fitting for /nl'"' in both cases. 

• We have shown that the systematic effect of stars on the DR9 
CMASS galaxy density field significantly affects the /nl'^' con- 
straints, but that the systematic bias has a different scale depen- 
dence than the (convolved) scale dependence of the effect of /nl^'- 
Thus, each effect is detectable in the data, and /nl^' constraints 
can be obtained even when allowing any amplitude of the system- 
atic correction (case iv), resulting in a 17% increase in the width of 
the 95% confidence interval. 

• We find that the data exhibit a marginal preference for a stronger 
systematic correction than provided by our fiducial weights for stel- 
lar density, as the minimum decreases by 2.9 when the correc- 
tion for stellar contamination is 45% stronger. We find no physical 
model that produces a similar change in the minimum but also 
find no additional systematic effect that can explain the need for 
such a systematic correction. 

We encourage all future studies to incorporate systematic un- 
certainties in a manner similar to that presented here in order to 
obtain /nl^' constraints that are robust to systematic uncertainties 
related to the treatment of the data. 
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APPENDIX A: WINDOW FUNCTION CONVOLUTION 

Here we outline the process of determining [kj], used in Eq. 

[8] to obtain the model we compare to the measured power spec- 
trum. 

The "spherically averaged power spectrum" within some k- 
volume Vfe can be modeled as the convolution between the true 
underlying power spectrum, and a window function, W 



Pm{k) = — / d'fc' 



(27r)3 



Pt{k' + e)W{e). 



(Al) 



We measure the power at particular k. Approximating th e fc-y olume 
as a thin shell, ^ d^k' = / dQkSnik - \k'\), Eq. 
rewritten 



Al 



can be 



Pmik) 



(27r)3 (47r) 



dO,k 



dn, / dr,Pt{k + e)W{e),(A2) 



where = |fc + e]. We change variables from (Q^, Qk) to (fl^ — 
Qk, f^e). Assuming isotropy, Pt{k + e) can be removed from the 
integral over fl^ and we obtain 



Pmik) 



(27r)3 (47r) 



dr, / d(n,~nk)Pt{k+e) / dn,W{e)(A3) 



We now define 
W,{e) = j dft,,W{e')So{r,, ~r,), 
which is independent of angle and obtain 



(A4) 



(27r)3 



dr,Ws{e) / d{n, - nk)Pt{k + e) (A5) 



and thus 



Pn.{k)^ I j^Ptik + e)Ws{^ 



(A6) 



Eq. |A6| is expressed in terms of the theoretical spherically- 
averaged power spectrum (e.g., output by CAMB), and the spheri- 
cally averaged window function. Thus we see that, because P(k) is 
spherically symmetric, the expected value of the measured power 
only depends on the spherically averaged window function, which 
we estimate from the random catalog and denote Pwin- Thus 



Pm{k) 



Pt{k + e)P„,„(6)e^(5(r,, - r,)dedcos(6l), (A7) 



where r^i = k + e and = ^Jk'^ ~\- — 2kecos{6). 

We measure the spherically averaged Pm and Pwin in discrete 
fc-bins. Further, determining | A7 1 f or every theoretical model tested 
would be computationally severe. To account for these facts, we fit 
a spline to Pwin and then determine the matrix via 



W[k,][k, 



{e)e^e{r„kj)dedcos{9) 



(A8) 



where now = + e2 - 2kiecos{e) and Q{r^, kj) is 1 if 
lies within the fe-bin kj and otherwise. Thus, we only need to 
calculate the input P(k) at each bin kj and use Eq.jsjto determine 
the model for our measurement at ki. Fig. [T] and the surrounding 
text demonstrates that this methodology works extremely well on 
mock DR9 CMASS samples. 



APPENDIX B: TESTS ON MOCKS 



The 600 mock DR9 CMASS samples created by [Manera et"ar] 
([20T2J allow us to test our methodology. We are using a covariance 
matrix constructed from the logarithms of the mock power spectra, 
P{k)mock- This procedure is tested by evaluating the skewness, G, 
of the distributions of P(k)mock and ln[P{k)mock] given by 



G(fc) = 



V60Q X 599 fX{ky~X{k) 



600 X 598 



E 



a(k) 



(Bl) 



where a is the standard deviation, X represents either P{k)mock 
or \n[P{k)mock\, and X represents the mean across all mocks. Fig. 
|B1 1 displays the skewness of the distributions of P{k)mock (red) 
and \n[P{k)mock] (black) for the fc-bins we consider in our anal- 
ysis. The expectation of the standard deviation (9.975 xlO~^) of 
the skewness of 600 values drawn from a Gaussian distribution is 
denoted by the black dashed line. For \n[P[k)mock], three of the 
twelve skewness are greater than the expected standard deviation 
and only the skewness at the lowest fc-bin is greater than two stan- 
dard deviations from zero. For P{k)mock, half of the skewness are 
greater than two standard deviations from zero. This suggests that 
testing ln[P(fc)] against models using the covariance matrix con- 
structed from \n[P{k)mock] is indeed a better choice than P{k) for 
the DR9 CMASS sample. 

We test our treatment of systematics for cases (i), (iii), and (iv) 
on the power spectra determined from the mocks (the mocks have 
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Figure Bl. The skewness, G, of the distributions of the 600 mock power 
spectra (red) and of the logarithm of the 600 mock power spectra (black). 
The dashed lines display the expectation of the standard deviation of the 
skewness of 600 values drawn from a Gaussian distribution (9.975X 10~^). 



Table Bl. Statistics for tests on mock DR9 CMASS catalogs: frac. 68; 95 
is the fraction of realizations that contain /J^l^' ~ ^ within the area of 
the PDF containing the highest 68.2%; 95% of the likelihood, J^^l^^ is 
the mean of the maximum likelihood value of /JJ'l^', cr is the standard 
deviation using /JJ'l^' = as the mean, and frac. w. 2 peaks is the fraction 
of PDFs that have two peaks. 



case 


frac. 68; 95 


Jlocal 


fj 


frac. w. 2 peaks 


i 


0.720±0.018; 0.960±0.008 


-14.2 


116 


0.090 


iii 


0.733±0.018;0.957±0.008 


-14.3 


116 


0.093 


iv 


0.688±0.018; 0.970±0.007 


-13.1 


120 


0.240 



within 1(7 of the expected value of 0.682, but the 95 fraction is 
2.9(7 from the expected value of 0.95. In order to bring the 95 frac- 
tion determined for case (iv) within Ict of 0.95, we must reduce the 
integration threshold to greatest 92.3% of the PDF. For both cases 
(i) and (iii), the 95 fractions are within 1.25(7 of the expected value, 
but the 68 fractions are greater than 2.1(7 from the expected value. 
In general, this suggests that the 95% confidence intervals we de- 
termine are more robust than the 68.2% confidence intervals, but 
that the width of the 95% confidence interval determined for case 
(iv) may be over-estimated. 



no systematic in their density field, so case (ii) is irrelevant). In 
general, we find that there is a degeneracy where in many case we 
find peaks in the /nl''' PDFs at both positive and negative values. 
This is due to the covariance: the best-fit value of /nl^' is most 
likely to occur when the model power spectrum has values in the 
two lowest-fc bins that are both either above or below the measured 
power spectrum. As summarized in Table |Bl| 9% of the mocks in 
cases (i) and (iii) have PDFs with two peaks and 24% of the PDFs 
for case (iv) have two peaks. We find that the mean /nl^' values 
(listed in Table |Bl| change by only 1.2 when we test the different 
cases, suggesting that our treatment of systematics does not bias 
the recovered /n^''' value. 

While consistent for each case, the mean /nl^' values are 
■~3(7 from the expected value of zero. This may be due to the 
fact that the distribution of \n[P(k)mock\ has a positive skewness 
(which implies median values smaller than the average) at the two 
lowest fc-bins. This potential bias is much smaller than the uncer- 
tainty expected on any individual realization. Most importantly, the 
bias is consistent in each case and therefore should not alter any 
conclusions we reach on the effect of observational systematics on 
the measurement of /n^'- 

We further test our methods by determining the percentage of 
mocks realizations that find the correct answer (/nl^' — 0) within 
68.2%/95% (which will denote as the "68/95 fraction") of the area 
of the individual PDFs, i.e., we integrate the regions of the PDF 
with the greatest probability until we include = and check 

if the integral is less than 0.682/0.95. We find that, given the prefer- 
ence for peaks at both positive and negative values, the best results 
are obtained when constructing the PDF of the absolute value of 
/nl''' when less than 90% of the PDF is on one side of zero (i.e., 
10% < P(/Nr' < 0) < 90%). The results of this test are in- 
cluded in Table |Bl| We find that the 68 and 95 fractions are higher 
than we expect for all cases, but that the behavior of case (iv) dif- 
fers from that of cases i and iii. For case (iv), the 68 fraction is 



APPENDIX C: CORRELATION FUNCTION 

Our study focuses on the results that can be obtained from mea- 
surements of the power-spectrum of DR9 CMASS galaxies. Since 
the correlation function and power spectrum contain the same in- 
formation, in principle, the same results should be achievable from 
measurements of the correlation function as our power spectrum re- 
sults. Here we present the analysis we performed on the correlation 
function. 

The spherically-averaged redshift space galaxy correlation 
function, ^o(s), is the Fourier transform of P°{k). Inspection of 
Eqs.|6] and |7] reveals, that, for a given background cosmology and 
redshift, four integrals over k are required to obtain ^o(s, 6, ./nl'^')- 
However, for the model £,o{s), we must also account for the non- 
linear smoothing of the baryon acoustic oscillation (BAG) feature, 
as the damping of high k BAO oscillations damps the (one) ^(s) 
BAO peak. We therefore use a damped matter power spectrum, 
PMd{k), obtained using a damping scale Sd (see, e.g., Seo &Eisen- 
|stein|2003||Crocce & Scoccimarro|2006[ l and: 



PMd{k) = PM{k)e 



(CI) 



Given PMd{k) and PM{k), we Fourier transform to determine 
the damped and undamped isotropic 3-dimensional real-space cor- 
relation functions ^Kn,d('") and ^un. We also include the effects 
of coupling between low and high k modes (see, e.g., |Crocce &] 
|Scoccimarro|2008[[Slnchez et al.|2008) via 

^(r)=e™,d(r)+A™4L^(r)eL(r), 



where is the derivative of (^un and 
1 



.(1) 



27r2 



PM{k)ji{kr)kdk. 



(C2) 



(C3) 



Given ^(r), one can calculate the correlation function terms analo- 
gous to the power spectrum terms in Eq.|7] 

We measure the correlation function, ^meas(s) for the ob- 
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Figure CI. The measured coirelation function of the DR9 CMASS sample, 
multiphed by (black) and the best-fit model, with f^ff^ = 218, shown 
by the solid red line. The black dashed line denotes ^(s) = and the dashed 
red line displays a model with /^l^' = 0. The baryon acoustic oscillation 
feature is seen at ~ 100/i~ ^Mpc. The significance of the apparent peak at 
~ 200?i"^Mpc is studied in Ross et al (2012). 
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Figure C2. The normalized (so that they integrate to 1) probability distri- 
butions for the local non-Gaussianity parameter /J^^^', obtained from the 
power- spectrum (P(fc); blue) and the correlation function (5(s); red), both 
using the fiducial systematic correction applied to the DR9 CMASS sam- 
ple. The distributions have similar maximum likelihood values for /J^l^'' 
but the ^(s) measurements allow a wider range of 



served DR9 CMASS sample using the standard [Landy & Szalay| 
( |1993^ estimator given by 

DD{a) - 2DR{s) 
RR{s) 



+ 1, 



(C4) 



where D represents the data sample (i.e., BOSS galaxies) and R 
represents the random sample (occupying the angular footprint and 
with the same redshift distribution as the data sample) and the pair- 
counts are normahzed to the total number. We use a hnear spaced 
binning in s of 7 h~^Mpc, and the measurements are the same as 
those presented in |Ross et al.| ( |2012) . 

By definition, the correlation function integrates to zero: 



^{s)ds = 



(C5) 



When estimating the correlation function, this condition holds, but 
we do not have an infinite volume. Instead the condition becomes 
a sum over the pairs from which ^meas is determined 

J:,RRis.)as^y. 



j:rr{s. 



= 0. 



(C6) 



Thus, for any particular model ^{s)mod, we can compute the devi- 
ation from this condition as 



j:RR{s)as). 



(C7) 



and thus the total model to compare to the measurement is 
^{s)m,od + Imod- For our fiducial model, and a bias of 1.9, it is 
Imod = 5.6 X 10-^ but for fj^l^' = 100, it is Imod = 1.3 X 10-^ 
(for reference, the standard deviation of the ^(s) computed from 
the mock catalogs is an order of magnitude larger, e.g., it is 0.002 
at s = 2Q0/i"^Mpc). 

Using the correlation function (as opposed to the power power 
spectrum) measurements presents additional challenges for the co- 
variance matrix that is applied. The amplitude of the measurements 
are generally expected to cross zero just beyond the BAO scale, and 
where they are close to zero we should not expect the fractional er- 
ror to be conserved for alternative models. We thus use ^(s) and 
its covariance, as determined from the 600 mocks, accepting that 
this may be inaccurate for large fj^f,^^, but that the results will still 
be useful for comparison. Further, the model ^(s) has two nuisance 
parameters, Sd and Amc, that the power spectrum does not have. 
For the results we present, we marginalize over these parameters, 
as in |Sanchez et al.| ( [2012| l, but we find less than a 10% variation in 
the width and maximum of likelihood distribution when we fix the 
damping scale to the value used in |Anderson et al.| ( |2012[ l. 

Fig.|cT]displays the measured for the DR9 CMASS sam- 
ple and the best-fit model ^o, applying case (i), for which we find 
f^l""^ = 218, consistent with the value of = 215 we 

found using the power spectrum. Fig. |C2| displays the probability 
distribution we find for /nl^' using the ^o(s) measurement with 
30 < s < 400/i^^Mpc (red) and the Po(fc) measurements at 
k < 0.05/iMpc~^ (both for the measurements using the fiducial 
case i). Each prefer positive /nl^' at greater than 68% confidence, 
but allow /nl"' < at greater than 5%. As expected, the two es- 
timators recover consistent results. However, the probability distri- 
bution is considerably wider for the correlation function measure- 
ments. This is likely due to the fact that the maximum effective 
scale probed is much larger for the power-spectrum measurements, 
for which the data compression at large scales is more effective. 
This result justifies our choice to focus on the power-spectrum re- 
sults. 



